Quark mass uncertainties revive KSVZ axion dark matter 
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The Kaplan-Manohar ambiguity in light quark masses allows for a larger uncertainty in the ratio of up to 
down quai'k masses than naive estimates from the chiral Lagrangian would indicate. We show that it allows 
for a relaxation of experimental bounds on the QCD axion, specifically KSVZ axions in the 2 — 3 /leV mass 
range composing 100% of the galactic dark matter halo can evade the experimental limits placed by the ADMX 
collaboration. 
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I. INTRODUCTION 

In the energy budget of the universe, it appears that approx- 
imately one quarter of the cosmos is composed of unidentified 
weakly interacting particles. The as-yet-unseen axion, origi- 
nally proposed to solve the strong CP problem, has couplings 
to normal matter suppressed by a large energy scale. As such 
it provides an ideal dark matter candidate. There are two pos- 
sible ranges of axion masses that would result in sufficient 
axionic matter density in the early universe to provide most, 
if not all, of the currently seen O dm- A mass scale of ~ jieV 
would provide cold dark matter through coherent oscillation 
of the axion field, while an eV mass would produce axions as 
thermal relics [1]. 

From measurements of the neutron dipole moment, it is 
known that the strong interaction preserves CP symmetry. 
This presents a conundrum for quantum field theory, as no 
gauge or global symmetry prohibits the CP violating d-tena 
in the Lagrangian: 



6'G''G". (1) 



This term is a total derivative, and therefore one is tempted to 
ignore it. However the topological structure of SU (3) allows 
Eq. (1) to contribute through instanton effects. Indeed, the ex- 
perimental fact that there is not a meson with mass < a/S^Itt 
which can be identified as the Nambu-Goldstone (NG) bo- 
son of a spontaneously broken U{1)a global symmetry of the 
strong interaction implies that there is no axial synmietry, and 
therefore we must take the presence of the 9 term seriously 
[2]. 

Experiments limit 9 ^ 10~^° [3] [4] [5], in contrast to naive 
expectations that 9 should be 0(1). The mystery deepens 
when one considers that the physical observable is not 6 itself, 
but 9 = 9 — arg det M where M is the quark mass matrix. It is 
difficult to see how 9, originating with the strong interaction, 
could cancel to such high precision with the phases of quark 
masses, originating in the electroweak sector. 

A favored approach to ehminate Eq. (1) is to seek a dy- 
namical method to set 9 to zero, by introducing a new de- 
gree of freedom: the axion [6] [7] [8] [9]. In brief, an axion 
model requires a new global Peccei-Quinn (PQ) symmetry, at 
least one scalar field, and strongly interacting fermions with 
PQ charges. The scalar field acquires a vev, breaking the 
PQ symmetry, and the resulting Nambu-Goldstone boson be- 
comes the axion. Through triangle anomalies, the axion mixes 



with mesons (resulting in a non-zero mass), and thus couples 
to photons, nucleons, and leptons. 

The axion mass and couplings are suppressed by powers of 
the axion decay constant /„. From astrophysical constraints 
/a > 4 X 10* GeV [10]. For fa ~ lO^^ GeV, the axion could 
provide cold dark matter with il ^ 0(1) with very weak cou- 
pling to the Standard Model particles. Many searches for ax- 
ions have been performed, so far all have retumed negative 
results [53]. 

We focus here on recent bounds placed on the axion to pho- 
ton coupling by the ADMX experiment [14], and the derived 
limits on the axion decay constant fa- The ADMX exper- 
iment is of particular interest as it has placed the most re- 
strictive bounds on the axion-photon coupling in the regime 
fa 10^^ GeV. Assuming that the dark matter in our galac- 
tic halo is comprised completely of axions and has a den- 
sity Pdm = 0.45 GeV/cm\ ADMX excludes the KSVZ 
axion model with E/N = [54] in the axion mass range 
1.9 — 3.3 /xeV. However, it was noted [17][18] that uncer- 
tainty in the masses of the light quarks due to the Kaplan- 
Manohar ambiguity [19] introduces larger than predicted un- 
certainties in the axion-photon coupling. We shall show that 
this allows the KSVZ model to evade the experimental bounds 
from ADMX. 

In section II, the full effects on axion mass and coupling 
formulae from the Kaplan-Manohar ambiguity are calculated 
for general axion models. Section HI then analyzes what re- 
gions of parameter space may be excluded by experiments in 
light of these ambiguities. We conclude in section IV. Details 
of the calculations as well as explicit formulae are provided in 
appendix A, while the analysis of two quark flavors is shown 
in appendix B. 



11. AXION THEORY 

We start with the most general possible axion model [55]; 
containing at least one complex scalar field a and with chiral 
femoions ■0 charged under a new global PQ symmetry. At least 
some of the tp must have non-trivial SU{3) quantum numbers. 

For a specific example, we consider the KSVZ model 
[15][16], in which the scalar field a has PQ charge +2, and the 
fermions are new SU{3) triplets. Vl has PQ charge +1, and 
ipR has charge —1. The assignment of PQ charges precludes 
explicit mass terms for the fermions, but allows Yukawa inter- 
actions ^Lcril^R. In this model, the SM fields are not charged 
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under PQ symmetry. Another popular model is the DFSZ ax- 
ion, proposed by Dine, Fischler, and Srednicki [21] and sep- 
arately by Zhitnitsky [22] . Here two Higgs doublets provide 
mass to the Standard Model particles, with another new scalar 
field to break the PQ symmetry. The SM fermions, as well as 
the three scalars, are all charged under a global PQ symmetry. 
As the SM fermions naturally fall into complete SU (5) multi- 
plets, this is a simple example of an axion model with particle 
content consistent with grand unification (GUT). 

We shall perform exphcit calculations in the KSVZ model, 
and note that the low energy axion-photon phenomenonology 
for DSFZ is identical with the parameter E/N (which we in- 
troduce shortly) set to 8/3. In the KSVZ model, the high en- 
ergy Lagrangian for fields with PQ charge is 

^ Kin. - Ki^Lcri^R + h.c.) - V{a) - ^OG'^G'^. (2) 

The scalar a acquires a vev w/\/2 through some dynamics 
of the scalar potential. The fermions get PQ-breaking mass 
terms, and the remaining massless mode ^ of the scalar a ap- 
pears as a phase to this term: 

H^LCT^R + h.c.) ^ ■^hviijLi^Re''^/'' + h.c.) (3) 

By judicious choice of a chiral rotation on ^, we can eUminate 
the angle in Eq. (2). To do so, we recall from the Fujikawa 
measure of the path integral that a rotation of e^"'^^ on a Dirac 
spinor contributes 



87r2 



(4) 



to the Lagrangian (here the are the generators of the 
representation of SU{3) to which ip belongs). Defining 
2 J2 Tiif^t^) = N6°''', where the sum runs over all if) with 
PQ charge (N = 1 for one Dirac S';7(3) triplet), and the ax- 
ion decay constant fa = v/N,-we see that a rotation of 



exp - 



(5) 



sends 9^9 + 4>/ fa in Eq. (2). It is important to note that, 
should any of the if) have C/(l) em charges, the transformation 
Eq. (5) will also add a FF term to the Lagrangian: 



C- 



E 



327r2 \N) fa 



-FF. 



(6) 



Here E = 2 ^ Q^, where again the sum runs over all PQ- 
charged Dirac fermions. Since i> is not a constant, this term is 
not a total derivative, and so cannot be ignored. 

As will be shown, at low energies (j) will get a vev, {(j)) = 
— faS, eliminating the constant 9 term in the Lagrangian. The 
axion then is the excitation of the field, a = <p — (0). At 
energies far below /„, we can integrate out the heavy degrees 
of freedom, leaving the effective Lagrangian 



FF. (7) 



To derive the axion coupling to photons, we eliminate the 
coupling of axions to gluons through rotation of the light 
quark fields. Later, as we move below the QCD scale, this 
will result in mixing between axions and NG mesons of the 
broken chiral SU{2) x SU{2). Exphcitly, we rotate 



exp 
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(8) 



(where q = u, d, s). The factor of three is the number of 
quark flavors being rotated, with an additional factor of two 
for left/right handedness. 

After the rotation Eq. (8), the quark-axion sector of the La- 
grangian Eq. (7) is 



+ (gLMe"/3^»te + h.c.) 
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(9) 
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here M = diag(m„, md, m,g) is the light quark mass matrix. 

We now consider the Lagrangian below the QCD scale, 
where the quarks have hadronized into mesons. In doing so 
the kinetic mixing terms vanish, as we have (by design) ro- 
tated all three light quarks by the same phase and there is no 
flavor singlet in the chiral Lagrangian. 

To determine the axion-photon coupling and the axion mass 
we use the chiral Lagrangian, which contains powers of the 
quark mass matrix M. To include the axion in this, we take 
our lead from Eq. (9) and simply rotate M by the phase 
^laj^Sa iu tijg Lagrangian. 

Let us now consider the form of the chiral Lagrangian, ig- 
noring for the moment the axion phase rotation. To leading 
order in M, the mass term is simply 



-Cmass = ^/x/2Tr(MS) 



h.c. 



(10) 



where E = cxp[2i7rT"//^] (a = 1, ... 8) is the meson field, 
T" are the generators of the adjoint representation of SU{2), 
jj, is an undetermined constant, and f^^ = 93 MeV. As is well 
known, one can expand out the exponential S field and, by 
taking the terms of 0(7r")^, find the masses of the NG me- 
son octet in terms of fi and the light quark masses (Eqs. (Al) 
and (A2)). Comparison with experiment then allows one to 
measure the ratio of masses, m^/md = z and niu/ms = w. 
The resulting value for z {z = 0.56) we shaU refer to as the 
Weinberg value for mu /md [23]. 

Adding in the axion field, the procedure is identical. One 
simply expands both S and e"*/^-'^" to second order in the ax- 
ion and meson fields and reads off their masses directly. As 
will be shown in more detail later in this paper, doing so re- 
sults in axion masses and couplings in terms of /x, z, w, fa, 
and E/N. Using the Weinberg value, one can then place lim- 
its on the scale /„ from null results of axion searches. 

It is also at this point that the (j) field (of which the axion is 
the excitation) develops a vev. Expanding S and considering 
the constant term, we see that the (j) potential is 
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+ h.c. 
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-M/^Tr(Af)cos 
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(11) 



This is minimized when cj) ~ {(p) ~ —Ofa, which justifies our 
expansion about this field value in parametrizing the axion a. 

Returning to the chiral Lagrangian, there is a significant 
ambiguity in the value of z that is not apparent from consid- 
eration of Eq. (10) [19]. A naive uncertainty on the Weinberg 
value of z would be 0.56 ± 0.05 [24]. However, the chiral La- 
grangian contains terms of higher order in M. Of particular 
interest is the Kaplan-Manohar term; 
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TrME + — Tr[det(M)M-iS^]^ + h.c. 

TO, 



This term approximately remaps 



(12) 



(13) 



(md and to, have similar transformations) where (3 = jiS/ms. 
By naive dimensional analysis, the Kaplan-Manohar term 
should have a coefficient jjp' f'^/K^, where — 'infr; is the 
strong scale [25]. Comparison with our choice of normaliza- 
tion gives 



(14) 



As /iTOs ~ 0.24 GeV2, 5 is - 0.2. 

As shown in appendix B, in the limit of two flavors, the 
dominant effect of the Kaplan-Manohar term is the remapping 
of Eq. (13), while it's inclusion in the three flavor case has 
more complicated results. This remapping means that the ratio 
z is essentially unconstrained by considerations of the meson 
masses alone. 

While the chiral Lagrangian alone cannot determine z, 
other methods are of more use, primarily lattice QCD and 
higher order chiral perturbation theory. The range on the ratio 
z adopted by the Particle Data Group is [26] [56] 



0.3 < z < 0.6. 



(15) 



The calculation of the range of z has considerable uncertainty 
associated with it, for example reference [27] quotes a value 
up to 0.8. To allow for this lack of a general consensus in 
the community, we therefore also include a more conservative 
estimate of the z range in our calculations: 



0.2 < z < 0.7. 



(16) 



This corresponds to a 20% contribution from next-to-leading 
order terms in the chiral Lagrangian [19], and is consistent 
with several measurements of the mass ratio [27] [28]. 

Once z is determined, the ratio w, parameters 5 and /iTO„ 
can be calculated as functions of z and the physical meson 
masses (Eq. (A6)) [19]. As detailed in appendix A, for sim- 
plicity we expand the relevant meson masses to leading order 
in w, exact solutions differ by less than 3% from this approxi- 
mation over the allowed range of z. Additionally, comparison 



with reference [19] shows variation of up to 10% in the solu- 
tion for It; as a function of z. This can be attributed to addi- 
tional higher order terms included in the analysis of reference 
[19]. The solution of 5 as a function of z is shown in Fig. 1. 
As expected, we find that 5 is a number of less than order one, 
and the Weinberg value of z corresponds to (5 = 0. 



Weinberg z value 




FIG. 1: Coefficient 5 of Kaplan-Manohar term in the chiral La- 
grangian as a function of z, the ratio of up to down quark masses. 
The Weinberg value for z — 0.56 corresponds to 5 = 0. The func- 
tional form of S is given in Eq. (A6). 

Applying the rotation M Me'"/^-''" in Eq. (12) to in- 
clude the axion, we expand and take the terms quadratic in the 
fields. From this, one obtains the axion/neutral meson mass 
matrix M'^ (Eq. (Al)). 

To complete the calculation of the axion-photon coupling, 
we recall that both the 7r° and rj couple to photons through 
triangle anomalies, and, from Eq. (Al), have mixings with the 
axion. The couplings to photons for the light mesons are 



^G^^^n^FF + ^G^^^rjFF : 
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32^ 



FF. (17) 



Denoting the explicit boson-photon-photon couplings as Ga-yy 
(see Eq. (9)), G„-,-y and G^-y^ for the axion, tt*' and 77 respec- 
tively, the Lagrangian above the mass of the rj includes the 
following terms: 
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(18) 



The heavy 77 and 7r° can be integrated out at low energies, 
and one is left with an effective low energy coupling, ga-y-y, to 
photons involving only axions: 



(?a77 — Ga-yy 
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.(19) 
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Here, we use the same normalization as reference [29] (and 
the ADMX collaboration), defining the coefficient of aFF in 
the Lagrangian as ^ga^-y- The full expression for the coupling 
is given in Eq. (A8). As a function of z, w, E/N, fa and 6, 
it can be rewritten as a function of only z, fa and E/N, using 
Eq. (A6). Evaluating Eq. (A8) when (5 = 0, we recover the 
well-known axion coupling at low energies (see e.g. [30]): 



9afj 1 



a 1 
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(20) 



Most experimental searches for axions make kinematic as- 
sumptions which place direct bounds on the axion mass nia 
rather than fa- Taking the eigenvalues of the matrix Eq. (Al), 
one finds a relation between the mass and coupling constant 
of the pion and those of the axions: 



'ml fa = rn^n'->fiF{z,w,S). 



(21) 



Here, m^o is the 7r°7r° entry in the neutral meson mass matrix 
given in Eq. (Al), and the full form of F is given in Eq. (A7). 
With S — 0, F agrees with the well-known result (again, see 
[30]): 



F\ 



s=o 



{1 + z){l + z + w)' 



(22) 



Expressing 6 and w as functions of z, we find that F does 
not have a large numerical difference from Eq. (22) over the 
allowed range of z. A useful rule of thumb is that, when S = 
0, 



ma 



6 /ieV 



10^2 GeV\ 

Ta ) 



(23) 



Thus, we come to the final form for the axion-photon cou- 
pling, expressed in terms of the axion mass, pion mass, pion 
decay constant, z, and fundamental constants (see Eqs. (A7) 
and (A8) for the exact expression): 



{1 + z){l + z + w) 



2tt fyrm^o 
' E 24 + z + w 
N 31 + z + w 



(24) 



0{wd) 



Systematic errors in the bracketed expression are introduced 
by two sources: the 3% error caused by using a leading-order 
expansion in w for meson masses when solving for w and 
S, and the 10% discrepancy between our function of w and 
that in reference [19] due to our omission of additional higher 
order terms beyond Eq. (12) (see appendix A). However, as 
shown in appendix B, the precise values of w and 5 are of 
lesser importance than that of z in Eq. (24), and as a conse- 
quence the systematic errors are numerically negligible in the 
final result. 

Surprisingly, for a wide range of z in Eq. (16), the exact 
form of ga-yj (in terms of the Kaplan-Manohar coefficient S) 
agrees very well with the expression with S set to zero. This is 
shown in Fig. 2. As explained in appendix B, this occurs be- 
cause S is associated with a factor of mj^, and in the two fla- 
vor chiral Lagrangian both terms in Eq. (12) transform identi- 
cally under the axion rotation Eq. (8). All this is to say that the 



main contribution of the Kaplan-Manohar ambiguity to axion 
phenomenonology is simply to allow z to vary widely; it's in- 
clusion was not necessary when deriving axion-pion relations. 
However, this realization was only apparent in hindsight. This 
observation justifies the discussion in reference [18]. 
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FIG. 2: Qajj/'ma and Qa-y-y /mali^Q ^ function of z for E/N ■- 
8/3. For the experimentally allowed range of z (0.3 — 0.6 or 0.2 
0.7) both calculations give very similar results 



An important point is the relative minus sign between the 
E/N contribution to ga-yf, which has its origin in the quantum 
numbers of the fermions charged under Peccei-Quinn, and 
the terms involving z, w, and 6, which arise from the axion- 
meson mixing. For appropriate values of E/N, near perfect 
cancelation can occur between the two contributions, strongly 
suppressing the coupling. Using the Weinberg value for z, 
one finds « -1.92 with an uncertainty of ±0.08, 

sparking much interest in axion models with E/N ~ 2. Us- 
ing the range of z from Eq. (15), this term (depending on z, 
w, 6), is -2.20 at z = 0.3 and -1.89 at z = 0.6. Should the 
more conservative estimate of Eq. (16) be taken, this term can 
vary from —2.33 to —1.81. 



What are possible values for E/Nl The original KSVZ 
model used new heavy uncharged 'quarks' as the PQ fermions 
■0. Hence, this model has E/N — 0. This is sometimes re- 
ferred to as 'the' KSVZ model, but it should be emphasized 
that the model would not suffer intrinsic flaws if the new par- 
ticles were charged. If we assume grand unification such that 
the PQ fermions are in a complete multiplet of SU{5) (which 
occurs in DFSZ axion models [21][22]) then E/N — 8/3, as 
can be quickly verified by explicit calculation of the 5 repre- 
sentation of SU{5). Other values of E/N are certainly possi- 
ble, requiring only novel quantum number assignment. 



At z = 0, -il±4^ - 

3 l-\-z+w 



= 8/3, and so the coupling to pho- 
tons would vanish for GUT axion models with a massless up 
quark. However, with z = 0, the up quark can undergo a chi- 
ral rotation as in Eq. (8) without introducing a phase to the 
quark mass matrix M. Thus, when z = 0, the axion is not 
necessary to solve the strong CP problem. 
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III. EXPERIMENTAL BOUNDS 

With the resuhs of the previous section, we can now exam- 
ine the recent bounds placed on the axion-photon coupHng. 
We start with general arguments limiting fa, then consider the 
ADMX results. In the construction of the Peccei-Quinn res- 
olution to the 9 problem, there were no a priori assumptions 
about the scale of fa- Indeed, the original formulation of the 
axion model identified the scalar a as the Higgs field, and thus 
fa as the electroweak vev [31] [32], though this possibility was 
soon ruled out [33] [34]. General astrophysical and cosmo- 
logical considerations greatly reduce the possible parameter 
space for fa and nia [57]. Note that none of these constraints 
have been reevaluated in light of the full uncertainty on z. A 
full consideration of the Kaplan-Manohar ambiguity's effect 
on axion physics is forthcoming [36]. 

Axions would be produced in the early universe, and so 
their properties must be such that they do not contribute a 
mass density ft a > 0{1). There are two separate situations 
that must be considered. If the universe never reached a tem- 
perature T > fa after inflation, then axions are produced via 
coherent oscillations due to some initial misalignment of the 
axion field away from the vacuum. In this scenario, axions 
were never in thermal equilibrium with the rest of the matter in 
the universe, and so started their existence as cold dark matter. 
To avoid overdosing the universe, fa ^ 10^^ GeV, translating 
to ma > 1 AieV [37][38][39][40][41]. Thus, nia - fieW is 
very appealing, as they provide a source of cold dark matter 
with the correct 0(1) density. 

Note that this assumes a 'natural' size for the initial mis- 
alignment angle after inflation. Various arguments have been 
presented to avoid this bound: anthropics may make it more 
likely for observers to exist in a universe with an 'unnaturally' 
small initial misalignment angle [42], some strong dynamics 
at high energies could force the initial misalignment to be 
small [43], or late entropy production from particle decays 
could sufficiently dilute the axions to avoid overclosure [44] . 

Alternatively, if the reheating temperature after inflation 
was greater than fa, then topological defects would form as 
the universe cooled. These cosmic strings would radiate ax- 
ions, again these relics cannot overdose the universe. Limits 
from this scenario imply fa ^ 10^^ GeV, or nia ^ 10~^ eV 
[45] [46]. Note that models with iV > 2 have an exact Zn 
symmetry that is spontaneously broken, and hence lead to cos- 
mologically unacceptable domain walls. 

In addition to these cosmological arguments, axion prop- 
erties can be bounded by several additional considerations. 
These fall into two broad categories: direct searches and lim- 
its on novel energy-loss mechanisms from stars. These con- 
strain the axion decay constant on the scale fa < 10^ GeV 
(ma > 0.01 eV), and so are not relevant for considerations 
of axionic cold dark matter. For a recent review of these con- 
straints, see reference [10]. We note that the CAST collab- 
oration has released new results [47] which supersede those 
quoted in reference [10]. 

With this in mind, the ADMX experiment is very interest- 
ing [ 14] [48] [50] . ADMX is a laboratory experiment searching 
for axions in the galactic dark matter halo through their cou- 



pling to photons. Axions interact with a strong magnetic field, 
allowing them to transition to photons with energy equal to the 
rest mass plus small kinetic corrections. The conversion 
occurs within a resonant cavity, which can be 'tuned' to look 
for photons of a particular wavelength [49], corresponding to 
an axion mass between 1.9 and 3.3 /leV for the ADMX ap- 
peratus [50]. 

ADMX used a local dark matter halo density of 
0.45 GeV/cm'^; assuming that this is composed entirely of ax- 
ions, a lack of detected signal can be interpreted as an upper 
limit on the axion-to-photon coupling ga^-y- Alternatively, one 
could place limits on the percentage of the halo which is made 
of axions. The excluded region is nearly constant if plotted 
as ga-fj/ma versus nia [14]. However, there is considerable 
structure on small scales, as seen in Fig. 3 [50]. Furthermore, 
more stringent bounds are available from ADMX's high reso- 
lution search over the narrow mass range 1.98 — 2.17 eV [48]. 
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FIG. 3: The experimentally excluded region at 90% confidence 
from the ADMX experiment, taken from [50]. Our estimate for 
the experimental upper bound for all values of nia, gl^^^y/m'i < 
1.44 X 10~^^ GeV~'^/eV^, corresponds to the coupling evaluated 
at the Weinberg value z — 0.56. The theoretically predicted values 
of {ga-y~i /maf for various values of z are indicated. The z value of 
0.733 corresponds to the smallest coupling for E/N . The values for 
z < 0.3 lie above the plot range. 

Using the Weinberg value for z and assigning E/N = 
for KSVZ and 8/3 for DSFZ, the ADMX collaboration com- 
pared the experimental upper limit to the predicted coupling 
for these two common axion models. In doing so, ADMX 
ruled out the KSVZ model as the sole source of dark matter in 
the galaxy (at 90% confidence) in the mass range over which 
they scanned (see Fig. 3). However, we note that this exclu- 
sion occurs on the margins of the experimental limits [58]. 

In addition, there are three additional considerations in play. 
First, z has a larger range of possible values than previously 
considered. Secondly, KSVZ models can have E/N ^ 0. 
Lastly, the precise value of the local galactic halo density is 
uncertain, and may not be composed of only axions. The 
ADMX limit can therefore be considered an exclusion of a 
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'benchmark' KSVZ model, with E/N = 0, together with 
the extra assumptions that z = 0.56 ± 0.05, and poM = 
0.45 GeV/cm3. 

In Fig. 4 the E/N dependence of {ga-^-y/nia)^ is plotted for 
the full range of allowed z. The coupling is divided by nia so 
that the experimental limit is independent of m^, and squared 
to agree with the ADMX convention. For E/N = 0, the 
allowed range of z permits only a small reduction in the axion- 
photon coupling. However, as is more clearly demonstrated in 
Fig. 3, for z — 0.6 this reduction is sufficient to allow a small 
window in {ga-yy/iria)^ which has not yet been ruled out by 
ADMX. 




FIG. 4: [ga-yy/rria)^ versus E/N for varying values of z. The 
ADMX exclusion line is the exclusion line as in Fig. 3, and corre- 
sponds to z = 0.56 and E/N — 0. The value of {gaji /nia)'^ as 
a function of E/N for the Weinberg value, z — 0.56, is indicated. 
The orange band corresponds to 0.3 < z < 0.6, while the wider 
blue band is the range from 0.2 < z < 0.7. The minimal coupling 
at E/N = 0, corresponding to z = 0.733, is also shown. 

An alternate visualization is shown in Fig. 5, which demon- 
strates the large parameter space of E/N versus z which has 
not been excluded by ADMX. As can be seen, even with the 
Weinberg value of z, a sizable region is still experimentally 
allowed. Additionally, when E/N — 0, values of the light 
quark mass ratio 0.56 < 2: < 0.8 allow g^^^ to escape the ex- 
perimental bound. Thus, even with the more aggressive bound 
on z from Eq. (15), it is premature for us to conclude that the 
KSVZ model in the 2 — 3 iieV range cannot be the dominant 
source of dark matter in the galaxy. 



IV. CONCLUSIONS 

The /ieV axion mass scale probed by ADMX is especially 
interesting, as such axions would provide an ideal candidate 
for dark matter. In order to compare experimental results 
with theory, an accurate knowledge of pion physics and quark 
masses is required. As we have demonstrated, the Kaplan- 
Manohar ambiguity in quark mass implies that the canonical 
Weinberg value of z, used in the derivation of axion-photon 
couplings, is very suspect. 
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FIG. 5: Contour plot of {gaj-y/rria)^ in units of GeV~'^/eV^ as a 
function of z and E/N. Darker regions indicate larger values. The 
line marked "ADMX" is the experimental exclusion estimated from 
90% CL exclusion in Fig. 3, taken from [50]. 



The ratio z is poorly constrained from QCD physics, sig- 
nificantly expanding the theoretical uncertainty on the pho- 
ton coupling for specific models. Allowing z to vary within 
the range 0.3 — 0.6 (0.2 — 0.7) allows the axion coupling to 
photons due to meson mixing to vary by 10% (15%). This 
dwarfs the systematic error in our calculations due to addi- 
tional higher order terms in the chiral Lagrangian. 

Due to this uncertainty, even with a fixed value of E/N — 
0, the KSVZ model can escape the ADMX experimental 
bound, as can be seen in Fig. 4. Similarly, outside the range 
E/N sa 2.1 ± 0.2, the possible variation in the value of z 
generally lowers the measured combination (5077/^0)^ by 
factors of a few. 

In the region E/N w 2.1 ± 0.2, complete cancelation be- 
tween the E/N and the axion-meson mixing term allows the 
magnitude of ga-y-y to evade experimental bounds by many or- 
ders of magnitude, even reaching zero in some cases. Perfect 
cancelation requires fine-tuning, but it is notable that near per- 
fect cancelation occurs for a considerable range of E/N and 
z values, as can be seen in Figs. 4 and 5. 

The DFSZ model is essentially unconstrained by the 
ADMX results. As can be clearly seen in Figs. 4 and 5, even 
with the Weinberg value for z, models with E/N = 8/3 fall 
well below the experimental limit. Fortunately for the discov- 
ery prospects, the GUT value of 8/3 can only be cancelled by 
z = 0. In addition to being outside the range we consider, this 
value for z would entirely negate the need for the PQ mecha- 
nism as a solution to the strong CP problem. Therefore, one 
should remain optimistic for the prospects for future bounds 
on the DFSZ model. 
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It should also be noted that the ADMX limit excludes only a 
benchmark KSVZ model with E/N = 0. Naively, one might 
assume that a model with no explicit coupling between the 
axions and photons due to charged PQ fermions would be the 
most conservative choice when placing experimental bounds. 
However, as has been demonstrated, the smallest couphngs to 
photons instead occur when E/N « 2, as a result of cancela- 
tions in opposite sign chiral rotations on the PQ fermions and 
standard model quarks. 

Should the ADMX collaboration increase their sensitivity 
by at least an order of magnitude as projected [51], then ax- 
ionic galactic dark matter models with either E/N = or 
E/N = 8/3 would be experimentally accessible for all possi- 
ble values of z. For both of these points, near-perfect cancel- 
lation cannot occur for the range of z. 

As stated in section III, there are numerous other axion 
bounds besides ADMX. The large uncertainty in z from the 
Kaplan-Manohar term should have similar effects on these 
limits as well. Preliminary work indicates that the bounds 
should in general be loosened by factors of 0(1). As many 
astrophysical constraints depend on couphng to nucleons and 
electrons, which generically lack the possibiUty of cancela- 
tion (as with the E/N term), one would not expect regions of 
parameter space where the bounds can be reduced by many 
orders of magnitude. ADMX, in presenting the strongest lim- 
its on axion phenomenonology, is also the experiment where 
the uncertainty in z is most critical. 



APPENDIX A: MASS MATRICES 

Expanding out Eq. (12) and taking the terms second order 
in the meson and axion fields, one can determine the NG bo- 
son mass matrix and the axion-7r°-?7 mass squared matrix Ad"^. 
The elements of this matrix are 



mlo 


= l^mu 


mlo^ 


= l^mu 


m\ 


= A*m„ 




= A*rn„ 




= l^mu 


ml 


= IJ^mu 



V3z 

{w + 4:Z + wz + w{l + 4:10 + z)6) 
3wz 

{l-z){l + 26)U 

3zfa 



(Al) 



{2z - w - wz + 2w(l -2w + z)5)f^ 

SV^WZfa 

{w + z + wz + Aw{l + w + z)5)fl 
^wzfl 

Note that these are the entries in the mass matrix, not the mass 
eigenvalues of the physical particles. In addition, the masses 
of the other NG bosons are 



{l + z){l + 6) 



2 2 {W + Z){l + W5) 

mrfo = m^o = M'Tiu (A2) 

wz 



(m|.±)QCD = l^ruu 



{l + w){z + w6) 



wz 



The subscript QCD refers to the mass squared ignoring elec- 
tromagnetic contributions. To leading order, the EM contribu- 
tions to the charged pion and kaon are the same [59]. 

To determine 5, w, and /im„ as functions of z, three phys- 
ical observables involving the meson masses and mixings 
Eqs. (Al) and (A2) are constructed. The physical 7r° mass 
can be extracted from the mixing of -k^ and 



I 2 \ 2 m^o^ 



{l + z){l + 5) w{z-lf{5-lf 



Az'^ 



(A3) 



+ 0{w^) 



The electromagnetic contributions to the physical tt* and 
mesons can be eliminated at leading order by taking the mass 
difference: 

(m^± - m^±)phys. = ^imu " ^ " ^) + ^(«') (^4) 

For the third observable, we expand the K'^ mass to leading 
order in w: 

mto = iimu ( - + - + s] +0(w). (A5) 

\w z J 

For calculational purposes, we truncated Eqs. (A3)-(A5) to the 
indicated order in w. However, using the exact expressions 
only introduces an error of < 3% in the numerical solutions 
for w, 5, and /xrn„. 

Combining Eqs. (A3), (A4), and (A5), we solve for ^m„, 
w, and 5 in terms of z and the observable masses of 7r°, K°, 
and the charged meson mass difference. Two sets of solutions 
are obtained, one gives negative values for w, and so is dis- 
carded. Using the known meson masses, the remaining set of 
solutions is 



S{z) 
w{z) 

fj,mu{z) 



1.78(-0.56 + z) 

-1.78 + z 
0.0279(z- 1.78)z 



^2-1 



(A6) 



(81.1 MeV) 



^{z-1.78)z 



1 



Fig. 1 shows (5 as a function of z. Comparison of our function 
w{z) with that given in reference [19] reveals ~ 10% dis- 
crepancy, originating in our neglect of additional higher order 
terms in the chiral Lagrangian. 

Diagonalizing the neutral meson mass squared matrix M'^ 
allows one to find the physical masses for the tt", tj, and axion 
a. Writing the axion mass in terms of the pion mass, we find 
the functional form of F{z, w, 5), introduced in Eq. (21). 

F{z,w,6) = ^[z+{w + z + wz)Sf + Azw'^S^^ 
/ {(1 + z)(l + S)[{z + z^ + wS){l + w + wS) 
+z{S + z6 - zw - wS + w'^d + w6^ + w^S"^)]} (A7) 
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Recall that w and S are functions of z. 

As in Eq. (18), diagonalizing the mass matrix also deter- 
mines the low-energy axion-photon coupling ga-y-y in terms 
of the explicit couphngs of a, 7r° and rj to photons, as given 
in Eq. (19). Using the expUcit mass squared terms from 
Eq. (Al), the full form of ga-y-y is found to be 



C = {z{4 + w + z) + ['w'^{l + z) + z{4 + z) + 



(A8) 



w{4 + z^)]S - 2w[w{z - 2) - 2z]S'^} 

/{{z + z^ + w5){l + w + w5) 

+z{5 + z5 — zw — w5 + vP5 + w5^ + w^5'^)} 



APPENDIX B: TWO FLAVOR ANALYSIS 

When working with only two quark flavors, including the 
axion proceeds identically up to Eq. (8). At this point, to 
avoid kinetic mixing terms in the chiral Lagrangian, the up 
and down quark fields are rotated by 



exp 



w 75 
fa 4 



(Bl) 



as there are only two quark flavors being rotated (four 
fields altogether, once handedness is included). Follow- 
ing the three-flavor analysis, we arrive at Eq. (12), with 



M = diag(mu . nid) and axions included by rotating M by 
exp(ia/2/a) (twice the rotation in Eq. (Bl)). 

At this point, it is clear why 5 enters into the axion phe- 
nomenonology only in the form wS. Naively, one would ex- 
pect ga-y-y to have terms of order S, wS, S"^, etc. The two 
flavor analysis corresponds formally to sending tz; — > 0, so 
the leading corrections would be of the form S. However, the 
leading order term in Eq. (12) transforms as e*"/^-'", while 
the {det{M)M-^)^ terms transforms as e-'^Z-^'^e^''/^^" = 
g-ja/2/„ Including the hermitian conjugates, we find that 
both the leading order piece and the Kaplan-Manohar correc- 
tion both transform identically under the axion rotation. 

Therefore, the Kaplan-Manohar term can be considered as 
a simple redefinition of m„ and rud as far as the axion is con- 
cerned. Thus, there is no 0{5) correction to the axion cou- 
plings. The effect to the higher order term can be completely 
parametrized in terms of z. 
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